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Depletion layer properties have been calculated for diffused junctions in 
silicon and germanium as a function of reverse voltage and of diffusion 
parameters for the gaussian and the complementary error function distri- 
butions. These results bridge the gap between the linearly graded behavior 
generally exhibited by such junctions at low voltage and the step behavior 
exhibited at high voltage. For total depletion layer thickness and capacitance, 
the transition from graded to step junction behavior extends over about one 
decade of voltage. For depletion layer thickness on a single side of the 
junction, it extends over several decades. Depletion layer thickness and peak 
electric field are presented graphically as a function of voltage for a variety 
of junction depths and impurity concentration functions. The ranges for 
which the step and graded junction approximations are valid are apparent 
from these charts. The results were obtained by an analytical integration of 
Poisson's equation, and a subsequent use of the IBM 704 for a numerical 
evaluation of the transcendental equations obtained. 

I. INTRODUCTION 

The dependence of depletion layer properties on voltage is important 
in the design of many semiconductor devices. Some of these properties, 
such as total depletion layer thickness, are accurately predicted by the 
commonly used step and graded junction approximations over fairly 
wide voltage ranges. However, the ranges in which the approximations 
are applicable have not been established previously. The present work 
establishes the regions in which these approximations are valid and 
supplies data for the entire voltage range of interest. 

Both the complementary error function distribution and the gaussian 
distribution have been treated. Total depletion layer thickness, peak 
electric field, capacitance per unit area and the fraction of the depletion 
layer on each side of the junction have been calculated for wide ranges 
of voltage, junction depth, background impurity concentration and 
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surface concentration. Results are given for silicon and germanium and 
can readily be extended to other materials. 

II. METHOD — COMPLEMENTARY ERROR FUNCTION DISTRIBUTION 

For the case of diffusion into a semiconductor with the assumption 
of a concentration-independent diffusion constant and a constant volume- 
concentration at the surface, the solution of the one-dimensional diffu- 
sion equation is the complementary error function. The solution may 
be written 

C(x) =C erfc^=, (1) 

where C is the impurity density at a distance x from the surface, C the 
surface concentration, D the temperature-dependent diffusion constant 
and t the time of diffusion.* 

For the purpose of these calculations we have taken the case of a 
donor diffusion into a p-type semiconductor. To treat the depletion 
layer formed at the junction, when reverse bias is applied, we have made 
the following customary assumptions: 

i. ionization of donors and acceptors is complete; 

ii. the depletion layer is a region completely free of carriers separated 
by a sharp boundary from an electrically neutral region. 

The net density of positive charge, p, lying within the depletion layer 
can be written 

p = q(C- C B ), (2) 

where q is the electronic charge and C B is the background impurity 
concentration. Hence Poisson's equation can be written in MKS units 
as follows: 






( c ° erfc ^br<4 (3) 



where ^ is electrostatic potential, k is the dielectric constant of the mate- 
rial and eo is the permittivity of free space. 

As Fig. 1 indicates, the electric field vanishes at the depletion layer 
boundaries. Thus, by integrating Poisson's equation, we can write 



erfc u = 1 — erf u; erf u = —7= / exp(— or 2 ) da. 
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expressions for the field distribution on the two sides of the junction: 
On the left, 



Ei(x) - Eyixj - 01) = E x {x) 



9C0 r ( 

K€ Jx,- ai \ 



, C B f x 



)* 



(4) 



and, on the right, 

Etixj + 02) - E»(x) = -EM 



qC 



n-i-^)-. 



K€o J \ Co 

where xj , the junction depth, is defined by 

Cnerfc— 7= - C B = 0. 



%Dt 



(5) 



(6) 





Fig. 1 — Diffusion profile and corresponding electric field distribution. 
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Requiring that the field be continuous at the junction enables us to 
arrive at the desired relation between the two components of depletion 
layer thickness, e^ and a 2 . That is, 

E 1 (x j ,^,a^ = E 2 (x j ,^,a 2 y (7) 

At this point specific values can be chosen for Xj , C b /Cq and ai so that 
a 2 can be calculated by a machine using iterative methods. 

Next, the voltage drops on the two sides of the junction are obtained 
by a second analytical integration of Poisson's equation: 

Vi= -[' E i (x)dx, (8) 

V 2 = -f ' ' E 2 (x)dx. (9) 

For the situation chosen, reverse bias implies a positive voltage on the 
left-hand side with respect to the right. Therefore, the positive x direc- 
tion is in the direction of decreasing voltage, and the calculated voltages 
are negative. 

The integrations in (4), (5), (8) and (9) are carried out in detail in 
Appendix A. 

Finally, junction capacitance is calculated from 

C = -^-. (10) 

(L\ T «2 



III. METHOD — GAUSSIAN DISTRIBUTION 

For the case of a diffusion having as a boundary condition that a 
fixed amount of the diffusant is deposited in the surface from which the 
diffusion proceeds, the solution of the one-dimensional diffusion equation 
is the gaussian function 

C = Co<T x2/(4D,) , (11) 

with the symbols being defined in the preceding section. 

Again we take the case in which the donors are diffused into a p-type 
semiconductor, and make the same assumptions as before. For this 
case, Poisson's equation can be written as follows: 



[-*2/(4D«) _ Cb 
Co J 



VV= £•«--«»« _£M (12) 

K€ 



DIFFUSED JUNCTION DEPLETION LAYER CALCULATIONS 393 

By integrating Poisson's equation as before we can write expressions 
for the field distribution on the two sides of the junction: 
On the left, 

EM - E l (x i - ai ) = E 1 (x) = ^ f \e-^ iDl) - C /\dx (13) 

K€ J Xj-ai L ^0 J 



and, on the right, 

Eiixj + 02) - Et(x) = -E 2 (x) 



KCq Jx [_ Co J 



(14) 



where Xj , the junction depth, is defined by 

[CiT'*' 14 * - C B ] = 0. (15) 

Requiring that the field be continuous at the junction enables us to 
arrive once again at the desired relation between the two components 
of depletion layer thickness, eti and a 2 . This equation, (7), is treated by 
machine as before. 

Next, the voltage drops on the two sides of the junction are obtained 
by the second analytical integration of Poisson's equation, as indicated 
in (8) for the left side and (9) for the right side. 

The integrations that yield expressions for E x {x), E 2 (x), V\ and V 2 
are completed in Appendix B. 

In this work, values of 16.00 and 12.00 were used for the dielectric 
constants of germanium and silicon respectively. A simple multiplica- 
tion enables one to determine field, voltage or capacitance for any 
modification of values of these constants, or for other materials. 

The IBM 704 was employed for solving the field continuity equation 
(7) and for performing all the numerical computations. 

IV. DATA PRESENTATION 

Fig. 2 is a diagram showing total depletion layer thickness versus 
voltage for three kinds of impurity distributions in silicon: step with 
large concentration on one side, linear graded, and complementary error 
function. Using gaussian results instead of complementary error func- 
tion results would give qualitatively the same picture. The relationship 
is shown for two junction depths: Xj = 2 X 10 -3 cm and Xj = 10" cm. 
The surface concentration is 10 20 atoms/cm 3 , and the background im- 
purity concentration is 10 16 atoms/cm 3 . The solid curves marked "erfc," 
which were obtained from the calculated data, clearly mark the transi- 
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Fig. 2 — Total depletion layer thickness vs. voltage for step, linear graded 
and complementary error function impurity distributions in silicon. 



tion region between the step and linear distributions, becoming asymp- 
totic to these two curves at the higher and lower voltages respectively. 
It should be noted that, since the step junction curve is a function of 
background resistivity only, it remains the same for both junction 
depths. 

The intersection of the linear graded and the step junction lines is 
the point where the accuracy of either approximation is poorest. In- 
creasingly larger errors occur, of course, if either approximation is used 
for voltages beyond its point of intersection with the other approxima- 
tion. The curves shown here constitute extreme examples. For Xj = 10~ 
cm the step junction approximation becomes excellent at about 4 volts, 
and hence the step approximation can be used over almost the entire 
voltage range of interest. In the case of Xj = 2 X 1CT 3 cm, however, 
breakdown occurs around 100 volts, and therefore the graded approxima- 
tion is good over almost the entire voltage range of interest. 

Fig. 3 presents total depletion layer thickness versus the ratio of volt- 
age to background impurity concentration. Fig. 4 presents the corre- 
sponding data for the gaussian distribution. These charts are available 
for a wide range of values of the ratio C B /Co ; a complete set will be fur- 
nished by the authors upon request. 

The chart in Fig. 3 is precise for C B /Co = 10 _a , but one may use the 
chart for a range from 3 X 10 -6 to 3 X 10 -0 . Maximum errors resulting 
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Fig. 3 — Chart for use in range 3 X 10 -6 to 3 X 10 -6 , erfc distribution. 



c io 3 - 




Fig. 4 — Chart for use in range 3 X 10~ 6 to 3 X 10 B , gaussian distribution. 
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from the use of Fig. 3 over this range are 5 per cent if one wishes to 
obtain total depletion layer thickness for a specific value of voltage. A 
maximum error of 10 per cent in voltage may result when starting from 
a given value of total depletion layer thickness. These maximum errors 
apply for the complete set of graphs. It can thus be seen that the use 
of a particular graph allows one to shift values of C B and Co as long as 
the ratio C B /Co falls within the specified range, and that, by using the 
complete set of graphs, one can find applicable data for any value of 
C B with its corresponding Co . 

As an example, if one has a ratio of C B /Co = 2 X 10 _ with C B = 
10 15 , one would select the chart applicable for the range 3 X 10 -6 to 
3 X 10 -5 (Fig. 3). One can then note that the voltage reading of one 
volt is obtained at V/C B = 10~ 15 . Then, for 10 volts, V/C B = 10~ 14 and 
Xj = 2 X 10 -3 , one obtains a total depletion layer thickness of 4.7 X 
10~ 4 cm. 

The plotted data yield parallel straight lines in the left-hand portion 
of each chart. These lines correspond to the graded junction approxima- 
tion, with each line representing a particular gradient. At high voltages, 
the lines then converge onto the single straight line that corresponds to 
the step junction approximation. Thus, these charts define clearly the 
applicability ranges of the two approximations insofar as total layer 
thickness and capacitance are concerned. Oblique axes have been used 
to spread the curves conveniently. 

Since V/C B is used as a parameter in the charts, avalanche breakdown 
data cannot be readily superimposed. The problem can be appreciated 
by noting that breakdown voltage depends on background doping in a 
step junction, but all step junction data have been collapsed onto a 
single line in these charts. 

Fig. 5 shows the variation of ai/a to tai with the voltage function for 
the complementary error function distribution. Fig. 6 gives the corre- 
sponding curves for the gaussian function. The upper boundary of each 
chart corresponds to the graded junction, and the lower boundary cor- 
responds to the step junction. Note that here the transition range covers 
many decades. 

A plot of the peak electric field divided by C B versus V/C B is given 
in Fig. 7. This chart was made for C B /C = 10 -6 , for the case of the 
complementary error function. The usefulness of the chart can be ex- 
tended to include a range of C B /C from 10~ 4 to 10 -8 , where the maximum 
error for a C B /Co other than 10~ is 10 per cent and, in the majority of 
cases, the error would be less than 5 per cent. For the gaussian distribu- 
tion, the usable range of this same graph is from C b /Cq = 10~ to C B / 
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v/Cb in volts-cm 3 for germanium 
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Fig. 5 — Chart for use in range 3 X 10~ 6 to 3 X 10~ B , erfc distribution. 
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Fig. 6 — Chart for use in range 3 X 10" 6 to 3 X 10 6 , gaussian distribution. 
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Fig. 7 — Peak electric field. 

C = 10 -6 , with Ch = 10 14 and with errors of the same magnitude as 
mentioned above. 



V. CONCLUSIONS 

For any given diffused junction there is a voltage below which the 
graded approximation is excellent and a voltage about a decade higher 
above which the step approximation is excellent for predicting total 
depletion layer thickness and capacitance. It can be seen from Fig. 2 
that rapidly increasing errors are introduced when either approximation 
is used beyond its applicable range. 

The approximations are less useful in predicting the division of deple- 
tion layer thickness between the two sides of the junction, for Figs. 5 
and 6 show that this transition region occurs over most of the voltage 
range of interest. The reason for this difference, of course, lies in the fol- 
lowing consideration: As voltage across a diffused junction is increased, 
ai becomes smaller than the half thickness predicted by the graded ap- 
proximation, and a 2 becomes larger. Thus, some compensation occurs 
when total thickness, a,i + a 2 , is computed. 

The results obtained here are immediately applicable to the calcula- 
tion of collector capacitance for junction transistors and junction capaci- 
tor design and the like. In field effect device design, where single-side 
depletion layer thicknesses are important, the results are particularly 
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helpful. Further, they should provide additional information for the 
detailed study of problems such as avalanche breakdown. 
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APPENDIX A 

As shown in Section II, we must solve Poisson's equation, which in 
this case can be written 

x 



VV = - 



«0 



' 



qC f 1 _Cb_ erf 



C Q 



r-Y 



(16) 



All the symbols are defined as before. From (4) we have, for the electric 
field distribution on the left side of the junction, 



and, from (5), for the right side, 

m\ dx ' 



»(»)-£•/* |Yi 



-k)- eri V4Di 



(17) 



(18) 



where Xj , a x and a 2 are defined in Fig. 1. 
Integrating (17) by parts, we obtain 

^-efi 1 -§)•- -TIB 



J A M e - xiniDt) ~\ 



C \ x i ~ a l 

1 _ cj (xj ~ ai) - (Xi ~ ai) erf VM 

4Dt -(.ij- a,)2/(4DO 



(19) 



Peak field will exist at the junction. Letting x = xj in this expression 
and noting that 

1 -&--*:* (20) 
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we can write for the peak field 

«%)-g{(%-*)(-^-- f^) 

_ i / 4.DI r -xj*K4Dt) _ --(*/- oi)V(4DO|l 



(21) 



For the range of arguments used in these calculations, the error func- 
tion is very close to unity, e.g., 0.99999980. Because the IBM 704 can 
handle at most eight decimal digits, it is advantageous to work with the 
error function complement. Making this substitution, the final expres- 
sion for peak field becomes 

El ( Xl ) = g {(„ - a,) (erfc ^= - erfc ^=') 

(22) 



_ ■ / 4:Dt r -1,-2/(41)0 _ -(*,— <u)*/(4i»)ll 



Similarly, (18) can be integrated to give an expression for the field 
distribution on the right-hand side of the junction and, specifically, the 
peak field: 

aw-g{(* + *)(-fc^--*.a^) 

/ 1 

_ , /4Dt r -XjtHiDt) _ -(*,-+o»)»/(4DO-jl 



The desired relation between a,\ and a 2 can be obtained by noting that 
field continuity requires 

E 1 (x j ) = E 2 ( Xj ). (24) 

Thus, we obtain an equation in a 2 , 

F(a 2 ) = 0, (25) 

which will be solved by trial and error. From (22) and (23) it is evident 

that 

F(«.) = «2 erfc ^= - (a, + a 2 ) erfc -^j=- 

+ j/^ a -*f*.W»> _ [ (ai - *,) erfc ^^ (26) 

- a, erfc ^|= + y/M e -<^«/<**>] . 
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Fig. S — Method of obtaining approximation a T r ■ 

For rapid convergence on a 2 in the machine solution of (25), Newton's 
approximation was used.* For a trial value, a T , which is a first approxi- 
mation to o-2 , we simply set a T = «i . By referring to Fig. 8 it can be 
seen that the next approximation, a TT , can be obtained from 

F(a T ) 

(27) 



&TT — (If 



(IF 
da 



1 ' 

Jar 



where 



(IF 

da 



1. - 



= crfc vmt ~ erfc vw • 



(28 



By substituting a TT back into (27) in place of a T (i.e., by repeated 
applications of Newton's approximation), rapid convergence results. 
A criterion of accuracy is the expression 



£,(.!•;) - EqJxj) 



(29) 



where the approximate value of a 2 is used in evaluating E 2 (xj). 

To obtain the total voltage across the junction (applied + built-in), 
(8) and (9) of the text must be integrated. Equation (8) yields 



ir ?£<> /[7i C *\ x * x * ,.r •'' 



V* xe 



v-f" 



-j-2/(4DO , 

e ax 



(30) 



This was suggested by Miss M. C. Gray. 
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where k is the second square bracket of (19). The last term of (24) can 
be written in erf form, yielding 

eH im- eH vm- (3l) 

Substituting (31) in (30) and converting to erfc form, we obtain 

xj - oA 

., (32) 
+ (xj + a0 J^ e -^-^i^ - x . A/^ e -^i^\ 

In a similar way, (9) yields 



r>--t[(™ + i-i){*^r^W^ 



'■-S[(«+*-*)\-^-<4?) 

+ (xj - as) Jn± e -^+°^i«™ _ Xi A/V± e - x i i,{iDt) ~\. 



(33) 
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We must again obtain a solution for Poisson's equation, which in the 
case of the gaussian distribution can be written 

aV = ^(V*' /(4i,I) -&V (34) 

«o \ Co/ 

where the symbols are defined as before. 

Integrations corresponding to those in Appendix A are carried out to 
obtain the electric field distribution on either side of the junction. These 
field functions are as follows: 

*«-S[**(-3r--^tc) 

c B 1 (35) 

- 7T ( x ~ x i + a l) 



a 

and 

7 



•w-eCvas^^s--^ 



- 7T (* - */ - fl 2) 

Co 
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The expressions for peak field are found immediately by substituting 
Xj for re in (35) and (36). Setting the resulting expressions equal to each 
other gives an equation that establishes field continuity at the junction 
[(7) of the text]. 

Integrating (35) again gives the voltage on the left-hand side: 

r.--g{vas%(-^--=fc§) 

+ 2Dt[e-^ liim) - e -^ )mwt) \ - ^V*''" 4 ™ J . 
Integrating (36) gives the voltage on the right hand side: 

+ 2Dt[e- Ii2KiDt) - «-<■*+■■>■'<*">] - \aie- x ' 2 "- iDl) \ . 



(38) 
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